Quantum Elastic Net and the Traveling Salesman Problem 
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Theory of computer calculations strongly depends on the nature of elements the computer is 
made of. Quantum interference allows to formulate the Shor factorization algorithm turned out to 
be more effective than any one written for classical computers. Similarly, quantum wave packet 
reduction allows to devise the Grover search algorithm which outperforms any classical one. In 
the present paper we argue that the quantum incoherent tunneling can be used for elaboration of 
new algorithms able to solve some NP-hard problems, such as the Traveling Salesman Problem, 
considered to be intractable in the classical theory of computer computations. 

PACS numbers: 03.67.Ac, 03.65.Xp, 03.65.Yz, 75.60.Ch 



I. INTRODUCTION 

Quantum parallel computations strongly differ from 
parallel calculi on the usual computers in a sense they 
use the same physical processor for all parallel opera- 
tions. For example, if one wants to solve some property 
recognition problem one can prepare an initial state in 
the form: 
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where x describes a possible assignment to binary vari- 
ables, Xi, of the problem, x — (xx, X2, x n ), X4 = 0; 1, 
and the last bit in the array vo is left to store a result of 
calculations, f(x), which is also presented in the binary 
form: 

m = 1, 

if the set (x\,X2, ...,x n ) has a required property, and 

m = 

otherwise. Calculations are performed during time t 
when the quantum computer evolves under unitary evo- 
lution operator, 

\v a )^\v t ) = U t \v a ) = ^=Y / \xJ(x)). (1) 

x 

In that way the results of calculations for all possible 
assignments of the variables X\, x%, X n turn out to be 
simultaneously evaluated and recorded on the same phys- 
ical equipment. 

It is clear that this method of computation gives a 
great economy of computer's time and memory. Unfor- 
tunately, the described straightforward strategy is ham- 
pered by the following essential reason. Quantum me- 
chanics says that if you try now to read out the results 
of calculations, you will succeed in observing only one of 
them, corresponding to some particular assignment of x, 



and all other results will be irrevocably lost. Further- 
more, the typical state of affairs with hard combinatorial 
problems is that results of calculations for most of assign- 
ments of x are equal to 0, and there are very few of them 
(and namely they are of interest!) which correspond to 
f(x) = 1. According to quantum mechanics, for the vec- 
tor \vt) in form {1}, any particular result of observation 
of f{x) will occur with probability P(x) = 1/2™. It is 
definitely forbidden to pick out at will any particular re- 
sult of calculations on the stage of reading out. You can 
amplify the contribution of a result you want to know in 
the total wave packet only on the stage of the dynamical 
evolution described, as in the above given example, by 
the unitary evolution operator Ut, or in the general case, 
by more general operator taking into account dissipation 
and/or measurement processes. In other case, the prop- 
erty recognition problem can not be efficiently resolved 
(is intractable) since the probability of observation of a 
necessary result is exponentially suppressed as P ~ 2~ n . 

Actually, there is no difference between calculations 
by quantum computers and the stochastic (Monte Carlo) 
calculations by classical computers, if one does not use 
some distinctive features of quantum systems to increase 
the probability of the result you want to obtain. For ex- 
ample, the quantum interference was used by P. Shor to 
determine a period of unknown function in his famous 
factorization algorithm 1 . This application of quantum 
mechanics is similar to usage of the interference for de- 
termination of a period of a crystal by Bragg's scattering 
method. The Grover search algorithm increases the prob- 
ability of the desired event with the help of the quantum 
wave packet reduction^. For the above considered exam- 
ple of property recognition problem, it is necessary to 
obtain after dynamical evolution a wave packet at which 
the components with f(x) = 1 are present with proba- 
bilities much bigger than the components with fix) = 0. 

In this paper we give some plausible arguments in favor 
of possibility of application of the quantum incoherent 
tunneling for solving some NP-hard combinatorial prob- 
lems. 
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II. ANALOG VERSUS DIGITAL 
COMPUTATIONS 

The main difference between analog and digital com- 
puters consists in a manner of representation of num- 
bers in them. In digital computers, the digital (usually 
binary) encoding of numbers is used, e.g., the value of 
x may be represented as 100101. In analog computers, 
different physical values, such as electrical currents, an- 
gles of rotation of gear wheels, etc, are representatives 
of mathematical values. Because of this the digital com- 
puters are much more compact and their size, L, scales 
with the quantity, n, of bits used for representation of 
the numbers as L ~ n. Enclosing such numbers in ana- 
log computers requires much more large physical storage 
device, L ~ e", and, therefore, corresponding resources, 
such as energy, forces etc necessary for operations with 
the numbers are also exponentially large in this case. 

One more difference following also from the manner of 
number representation concerns the property of univer- 
sality: digital computers are ordinarily universal, analog 
computers are designed for solving some special problems 
only. For each particular digital computer, universality 
should be proved, i.e. every arithmetic, logical, and other 
operations should be constructed of the basic operations 
with digits. 

Usually merely computations that do not use resources 
that grow exponentially are of interest (so-called efficient 
computations) . For them the Strong Church 's Thesis was 
formulated^: Any finite analog computer can be simu- 
lated efficiently by a digital computer, in a sense that 
the time required by the digital computer to simulate the 
analog computer is bounded by a polynomial function of 
the resources used by the analog computer. This Thesis 
resulted from a great experience gained during elabora- 
tion of classical computers and has a clear sense: every 
classical physical process can be efficiently simulated with 
the help of the digital computers. 

Due to a great progress in development of digital com- 
puter techniques, analog computers are used now ex- 
tremely seldom. Nevertheless, they survived in a vir- 
tual form inside digital computers as analog algorithms 
and turned into powerful heuristic methods for solution 
of optimization problems. The analog algorithms inher- 
ited from their natural ancestors the correspondence of 
computational schemes to some real physical processes 
in nature. Thus, the steepest descent method is based 
on gradient equations often used for far-from-equilibrium 
physical system description^. It is used to find local min- 
ima. Simulated annealing describes the opposite situa- 
tion when a physical system is very close to the equilib- 
rium state at each moment of its slow dynamical evo- 
lution. It is an example of analog algorithms for mini- 
mization of a function called energy which imitates the 
process of physical annealing known in practice of crys- 
tal growing. It is experimentally discovered and theoret- 
ically understood that molten matter, to be subjected to 
slowly cooling, transforms into a crystalline state corre- 



sponding to the minimum of its free energy. Simulated 
annealing was the only method (in addition, of course, 
to the honest looking over all variants) which gave hope 
to find the global minimum of combinatorial problems. 
However, rigorous consideration showed that in fact the 
global minimum can never be reached^, because simu- 
lated annealing with cooling schedule 

Ti, T2, T n — > 0, 

requires an exponential computation time: 

t n ~ exp(A/T„), A — const. 

Main computation time losses take place at the lowest 
temperature reached during calculations and the mini- 
mum becomes more and more inaccessible at each step 
of cooling. Even a very simple minimization problem cor- 
responding to a biased double well energy function shown 
in Fig.l requires the cooling schedule 

T n ~ D I In n 

and an exponential computational time, t ~ exprt, and, 
therefore, is intractable in the frame of this algorithm. 
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FIG. 1: Minimization problem corresponding to a biased 
double well energy function requires the cooling schedule 
T n ~ D/lnn and an exponential computational time in the 
frame of the simulated annealing algorithm. 

It appears at first sight, one can somehow accelerate 
cooling if does not require reaching the energy minimum 
on certainty. Then, after several reiteration, the position 
of the global minimum could be selected among all lo- 
cal minima found. But this strategy is, of course, wrong 
because no predictions about possibility to recognize the 
global minimum could be done if the true cooling sched- 
ule is broken. 

Description of quantum systems often requires consid- 
eration of exponentially large number of different vari- 
ants, e.g. a numerical realization of the path integral 
approach, and in particular for this purpose R. Feyn- 
man suggested quantum computers would be especially 
effective 5 . For quantum computations no statements 
similar to the Strong Church's Thesis were formulated so 
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FIG. 2: Tunneling with dissipation allows particle to pene- 
trate through the potential barrier in a finite space of time 
and then to stay at the point of the global minimum. 

far and, therefore, both analog and digital approaches are 
still of equal interest. Below, we discuss an approaches to 
solution of the Traveling Salesman Problem by a quan- 
tum analog computing machine. 

III. INCOHERENT TUNNELING 

In this section plausible reasoning are given in favor of 
incoherent tunneling as an effective remedy for solution of 
minimization problems on the stage of low temperatures 
when the simulated annealing algorithm is inoperative. 
To this goal let us consider again a biased double well 
potential shown in Fig. 2 and a particle localized initially 
in the left local minimum at point q^. Due to quantum 
tunneling effect, the particle can penetrate the potential 
barrier and fall into the global minimum at q±. When 
the motion of the particle is accompanied by dissipation, 
there is a chance that the particle will stay at the global 
minimum and, hence, an optimization problem will be 
resolved. 

A probable picture sketched above may be rigorously 
grounded. In the case when thermal energy, fcgT, is much 
less than energy, fiui, of classical oscillations at the po- 
sitions of minima and the latter is, in its turn, much 
than the height of the potential barrier, Vban between 
the wells, 

ksT <C flUJ <C Vbar, 

the process of tunneling with dissipation can be described 
on the base of multi-instanton calculations in the frame- 
work of the imaginary time functional integral approach 7 . 
For temperature close to zero, the tunneling rate from go 
to qi is given by an expression: 

. 7tA 2 1 /CT\2«-1 

2 uj 1 (2a) \ujJ 
while the inverse transitions are suppressed, 
7(9i -> Qo) = 0. 



Here a — i]d 2 /2irh is dimensionless damping coefficient, 
d = qi — <7 is the tunneling length, 77 describes the force 
of friction, F = —r/dq/dt, and A is the tunnel matrix of 
undamped and unbiased system. The time of transition 
to the minimum, r, is finite now: 

h 

T ~ — . 

7 

Thus, at least for the double well, tunneling with dissipa- 
tion, or incoherent tunneling, can be effectively applied 
for solving the minimization problem in the vicinity of 
zero temperatures. 




FIG. 3: Particle moves down the tilted periodic potential hill 
with a constant speed. 

In the case of a tilted periodic potential, displayed in 
Fig. 3, it is also possible to fulfil calculations explicitly and 
show that the time of transition from one local minimum 
to another is finite. An average position of the particle 
may be expressed as£ 

(q(t)) =vt, 
where v is mean drift velocity, 

v = — tanh(/3ficr/2), /3 = l/k B T. 

7 

Therefore, the time of minimization is proportional now 
to the sum of the tunneling time through all potential 
barriers on the way to the global minimum. It is natu- 
rally to suggest that a similar picture remains in the gen- 
eral case when the slope of the "hill" is variable. Thus, 
the main question left to be discussed is how many local 
minima actually are on the way to the global minimum. 
In the case of the exponential number of them the prob- 
lem would remain open. 

The practice of solving optimization problems shows 
that usually there are much less of energy local minima 
than all possible (non minimal) values of energy. For ex- 
ample, number of local minima of the energy function was 
thoroughly investigated for the Hopfield Network? which 
was suggested as a mathematical model of the associa- 
tive memory in human brain. Here each local minimum 
H(x) is used to retain an information about x (called a 
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pattern) since for a small deviation from x, x — > x + Sx, 
the system returns back into x if it is navigated with the 
help of the gradient equations. It was shown that for 
random created patterns the number of local minima M 
can be estimated as 

M = 0.138 n, 

where n is a number of artificial neuron cells^. Thus 
M is only a linear function of the size, n, of the input. 
A similar result was obtained for an associative memory 
based on Q-state Potts-glass with biased patterns^. 

Further increase of the number of stored pattern up to 
a = 1 converts the Hopficld Network into a spin-glass 11 . 
The problem of finding the ground state in a spin-glass 
was studied ini^. This problem was shown to belong to 
the class of NP-hard problems both for three-dimensional 
case, and for two-dimensional lattice within a magnetic 
field. Infinite-ranged models of spin-glasses were consid- 
ered inii Numerical experiments have shown the total 
number of local minima increases as some small power of 
n, rather than exp(n). 

A drastic decrease of the number of local minima with 
different energies in the vicinity of the global one, for any 
system in the thermodynamic limit, follows also from the 
Nernst theorem for entropy, 

lim S(T) = 0, 

T->0 

if there are no gaps of function S(T) near T = 0^. 

It is a commonly accepted method to demonstrate an 
efficiency of a new quantum algorithm with the help of 
a classical computer simulation of corresponding quan- 
tum calculations. But in our case there is no need to do 
so, because this has been already done. Actually, some 
kind of incoherent tunneling was used efficiently many 
times as an optimizing procedure during simulated an- 
nealing. For instance, not only neighboring (in the 
Hamming sense) configuration of spin were checked dur- 
ing minimization, but also states which were 2, 3 and 4 
steps away. When energy of a trial state was found to 
be less than energy of the current state the system was 
transmitted to the trial state. This corresponds to the 
under-barrier transitions through the potential barriers 
with 1, 2 and 3 Hamming's steps of width. Of course, 
such a strategy is hampered for classical computers for 
more wider potentials barriers by a huge number of pos- 
sible trial states. Quantum incoherent tunneling turns 
out to be much more effective because it allows to run 
over all local energy minima only, without examination 
of all possible values of energy for all possible trial states. 

IV. ELASTIC NET APPROACH TO THE 
TRAVELING SALESMAN PROBLEM AND ITS 
QUANTIZATION 

The Traveling Salesman Problem (TSM-problem) is 
formulated as follows: given positions of cities, what is 



the shortest tour in which each city is visited once? It is 
evident that the total number of possible tours increases 
exponentially with the increase of number, n, of cities: 

N* ~ e Cn 

1 v tours ° ) 

and their sequential consideration requires computing 
time which increases faster than any power of n. There- 
fore problem is considered as intractable in the classical 
theory of computation. In this section we show how the 
problem can be solved, at least in principle, using the 
quantum calculation technique. We consider the ana- 
log approach, because practical construction of quantum 
digital computers is still confronted with serious difficul- 
ties. Therefore, it is quite possible that quantum analog 
computers will win the digital ones, despite the opposite 
situation in the field of the classical computations. 

Quantum analog computers for solving TSM-problem 
can be constructed on the base of the Elastic Net analog 
algorithm elaborated for the usual digital computers^. 
The algorithm is grounded on a discrete form of the gra- 
dient equation, 

djjj __ _ dF 
dt di/j ' 

with free energy, F: 

F=-ak 2 ^Y / e-^ l/2k2 +0kJ2\mi-m 2 - 

i j 3 

For At = 1 it is possible to write: 

This means that the algorithm directs the system to a 
local minimum of F in accordance with the steepest de- 
scent method. The explicit form of F was devised ir>i^ 
using a very clear physical picture. Positions of cities are 
described by Xi, points with coordinates lie on an elas- 
tic string. Each point moves under the influence of two 
types of force. The first moves it towards those cities to 
which it is nearest; the second pulls it towards its neigh- 
bors on the string, acting to minimize the total string 
length. 

The authors described an iteration procedure which 
consists in a gradual decrease of the value of a parameter 
K describing a force and a range of interaction between 
x*i and yj in such a way that at K approaching to zero the 
range is transformed from a big value to very small one 
and the force is strongly increased. After applying this 
iteration procedure, an initial string in a form of a small 
circle placed nearby the cities' centre of mass is converted 
into a salesman tour of a very high quality. Comparative 
analysis fulfilled in^ has shown that the performance 
of the Elastic Net algorithm is even of higher quality 
than that of Simulated Annealing. Therefore, one can 
use any of them, or some other described ini£, to find a 
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FIG. 4: The process of solution of TSM-problem by the do- 
main wall analog computer: an input, i, corresponding to a 
local minimum spontaneously decays into the final state, /, 
describing the global minimum. 

preliminary solution to TSM-problem as an input for a 
further quantum computation. 

Elastic string can be created as a physical computation 
device in the form of a mesoscopic domain wall in an an- 
tiferromagnetic thin film at low temperature. Tunneling 
of the domain walls in ferromagnetic and antiferromag- 
netic insulators through potential barriers (see, e.g.^) 
is described in a similar way as incoherent tunneling of 
particles, considered in previous section. The role of the 
barriers can play defects or other singular points in the 
lattice. The domain wall has a variety of dissipative cou- 
plings to the magnons, photons, impurities, defects, and 
phonons. Nevertheless, it was shown that large domain 
walls, containing up to 10 10 spins can behave as quan- 
tum objects at low temperatures^. Both theory and 
experiment revealed a finite tunneling time of domain 
walls through potential barriers^. This gives an oppor- 
tunity of devising an analog quantum computer using 



the incoherent tunneling effect. The process of solution 
of TSM-problem by this computer is shown schematically 
in Fig. 4. Firstly, one creates an input, i, corresponding to 
some local minimum found with the help of the classical 
computer, and then, after a time, obtains a final state, 
/, as the solution. Points in the picture denote "cities", 
which are some kind of defects in the lattice that pin the 
domain wall to necessary locations. They may be some 
implanted atoms with high value of spins, or strongly 
magnetized atoms. Energy of the domain wall per unit of 
length (and, therefore, tension of the string) can be reg- 
ulated by the material properties, or by a change of the 
mutual orientation of directions of spins on either side of 
the domain wall. It is clear that a value of tension should 
be small enough to prevent the string from a temptation 
to pass round some cities to minimize its length. Such 
a slack string limit corresponds to the strong shot-range 
interaction between the string and the cities considered 
as the last step of the iteration procedure in the classi- 
cal Elastic Net analog algorithm^. Therefore, quantum 
computer solution of the problem can be interpreted as a 
natural continuation of the classical algorithm after tak- 
ing into account the tunneling processes. It is an easy 
task to examine that at least for a square lattice with the 
four nearest spin-spin interactions tension of the string 
is independent on its orientation. Thus, setting aside 
purely engineering problems, one may conclude that the 
Traveling Salesman Problem seems to be solvable by the 
quantum analog computer. 
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